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ABSTRACT 

There is growing evidence that the majority of the energy density of the universe 
is not baryonic or dark matter, rather it resides in an exotic component with nega- 
tive pressure. The nature of this 'quintessence' influences our view of the universe, 
modifying angular diameter and luminosity distances. Here, we examine the influence 
of a quintessence component upon gravitational lens time delays. As well as a static 
quintessence component, an evolving equation of state is also considered. It is found 
that the equation of state of the quintessence component and its evolution influence 
the value of the Hubble's constant derived from gravitational lenses. However, the 
differences between evolving and non-evolving cosmologies are relatively small. We 
undertake a suite of Monte Carlo simulations to examine the potential constraints 
that can be placed on the universal equation of state from the monitoring of gravita- 
tional lens system, and demonstrate that at least an order of magnitude more lenses 
than currently known will have to be discovered and analysed to accurately probe any 
quintessence component. 
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1 INTRODUCTION 

The searches for supernovae at cosmological distances have 
proved very successful, providing evidence that, while topo- 
logically flat, the majority of energy in the Universe is in the 
form of an exotic component with negative pressure (Riess 
et al. 1999; Perlmutter et al. 1999). The recent identification 
of a supernova at z = 1.7 (Riess et al. 2001) has provided 
further weight to these claims (Turner & Riess 2001), which 
suggest that this component may differ from the classical 
cosmological constant A. Termed 'quintessence', or more col- 
loquially 'dark energy', this has an equation of state of the 
form P — wp, where P is the pressure and p the density. 
A w < — | opposes the action of gravity and drives the 
cosmological expansion to accelerate. Linder (1988a; 1988b) 
has examined the physical nature of various quintessence 
components; with w = equating to non-relativistic mat- 
ter (dust), w = | being radiation and w = — 1, a classical 
cosmological constant. More exotic components are; mass- 
less scalar fields w = 1, cosmic string networks w = — -|, 
and two-dimensional topological defects w = — |. As well as 

1 It should be noted, however, that the influence of gravitational 
lensing on SN1997fT needs to be fully addressed before its true 
cosmological significance can be addressed (Lewis & Ibata 2001; 
Moertsell, Gunnarsson & Goobar 2001) 



the supernova programs, other approaches, such as gravita- 
tional lensing statistics (Cooray & Huterer 1999), geometri- 
cal probes of the Ly a forest (Hui, Stebbins & Buries 1999) 
and galaxy distributions (Yamamoto & Nishioka 2001), and 
classical angular-size redshift tests (Lima & Alcaniz 2001), 
will provide complementary probes of the universal equation 
of state. 

The value of the quintessence component, w, influences 
our view of the universe, modifying the various distances 
used in mapping the cosmos. This paper concerns itself with 
the influence of w on angular diameter distances, especially 
in relation to the determination of the Hubble's constant 
from the measurement of time delays in gravitational lens 
systems. Unlike local determinations of Hubble's constant 
(e.g. Freedman et al. 2001), the cosmological nature of grav- 
itational lenses means that they are more sensitive to the 
underlying cosmological parameters. Section ^| briefly cov- 
ers the basic formulae for generalized angular diameter dis- 
tances in quintessence cosmologies, while in Section ^ we 
consider the influence of w on the determination of H from 
lensed systems. Section ^ extends this analysis to simple 
models of an evolving quintessence component. In Section ^ 
a series of Monte Carlo simulations are undertaken to es- 
timate the efficacy of this approach in probing the cosmo- 
logical equation of state, while in Section ^ we speculate 
on the possability that current observations of gravitational 
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Figure 1. The relative time delays for several flat cosmologies as a function of the equation of state, w. Each curve is is for the lowest 
value of the time delay (the redshift of the lens-source pair is given in brackets). Each panel presents a different combination of Q m and 
Q m ; note that w = corresponds to a universe composed entirely of matter in all cases. 



lens systems may suggest that w < — 1. The conclusions of 
this study are presented in Section |^ 



2 GENERALIZED ANGULAR DIAMETER 
DISTANCES 

While there has been a resurgence in quintessence cosmol- 
ogy, the generalized cosmological equations for such uni- 
verses were presented more than a decade ago by Linder 
(1988a;1988b), including a generalized form of the Dyer- 
Roeder equation for the evolution of a bundle of rays travel- 
ing from a distant source (Dyer & Roeder 1973). Expressing 
the angular diameter distance as D = (c/H )r, the general- 
ized beam equation is given by 



r + 



3 + q(z) 

1 + 2 
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where Q w (z) is the density, in units of the critical density 
p c (z), of a contributor to the total energy-density of the 
universe with an equation of state w = P/p, and where P is 
its pressure and p is its density. This is given by 



Q. w {z) = 



Pc(z) 



Ho 

H(z) 



(2) 



where f2m(0) is the contribution of this component to the 
present energy-density budget at the present epoch, and the 
critical density is given by p c (z) = 3H(z) 2 /8ttG. Here, H(z) 
is a generalized form of the Hubble parameter and is given 
by; 



H(z) = H 



K(l + z) 



(3) 



where K = f2 — 1, and fl = Q(0) — ^2 ^tu(O), is related 
to the overall curvature of the universe. q(z) is a generalized 
form of the deceleration parameter and is given by 
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Finally, Equation [t] also contains the parameter a w (z) which 
represents how much of the fluid lies in the beam and influ- 
ences the the evolution of a ray bundle. For a universe con- 
taining matter, the solution to Equation |l] with 0:0(0) = 1 
represents the classic Dyer-Roeder 'filled beam' distance, 
while ota(z) = is the 'empty beam' distance (Dyer & 
Roeder 1973). 

When solving Equation the boundary conditions need 
to be defined. These are 



dr(z , z) 



dz 



= (1 + Zo) 



r(z ,z ) = 
H 



H{z) 



(5) 
(6) 



Equation [j] was integrated using a Runge-Kutta scheme (the 
rksuite package from www.netlib.org) and compared to 
both analytic results and the minimum angular extent red- 
shifts in quintessence cosmologies as tabulated in Lima & 
Alcaniz (2000); excellent agreement was found. Throughout 
this work, filled-beam distances a w (z) = 1 are employed. 



3 GRAVITATIONAL LENS TIME DELAYS 

Refsdal (1964) was the first to note that cosmological pa- 
rameters could be determined from the measurement of a 
time delay between the relative paths taken by light though 
a gravitational lens system. Since the discovery of multiply- 
imaged quasars, this has become the goal of a number of 
monitoring campaigns (e.g. Cohen et al. 2000; Oscoz et al. 
2001; Patnaik & Narashima 2001), although these analyses 
are frustrated by degeneracies in the derived mass models. 
The cosmological model simply enters the determination of 
Hubble's constant; 



H At oc 



TolTos 
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where m are the normalized angular diameter distances be- 
tween and observer (o), lens (I) and source (s), and At is 
the measured time delay between an image pair. 

Giovi & Amendola (2001) examined the influence on 
a quintessence component on gravitational lens time delays 
and the determination of Hubble's constant. Their analy- 
sis, however, was mainly concerned with the influence of the 
clumping of material and the dependence on H of whether 
distances are empty-beam or full-beam. Here, a different 
approach is considered; Solving Equation |l|, Equation [?] is 
evaluated for a range of quintessence components. For the 
study, the redshifts of the lens and source in seven gravita- 
tional lens systems that are favourable for time delay mea- 
sures were considered [Table 2 in Giovi & Amendola (2001), 
Q0957+561 and B1608+656], plus two fiducial redshift pairs 
of (zi = 1.5, z a = 2.5) and (zi = 1.0, z s = 2.0). While there 
are currently no lensed systems with established time de- 
lays at these particular redshifts, there are several potential 
systems; e.g. the quadruple lens H1413+117 at a redshift of 
2.55, with a lens redshift, established from prominent ab- 
sorption features, at ~ 1.5 and HE 11 04- 1805 at a redshift of 
2.31 and an estimated lens redshift of z ~ 1. ft is assumed 
throughout that the universe is flat, Q m + Q m = 1. 
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Table 1. The minima for the curves in Figure [l] for the various 
cosmologies under consideration. 

Figure ^ presents the results of this analysis; six pan- 
els are presented, each for a different combination of Q m 
and fl w . A greyscale-coded line for each redshift pair is pre- 
sented. The abscissa presents the equation of state param- 
eter, w, while the ordinate presents the relative time-delay; 
this represents the change in the time delay for a fixed Hub- 
ble's constant. Conversely, this is the relative value of H for 
a system with a fixed time delay. Each curve is normalized 
to the minimum value of the time delay. The relative time 
delay depends quite strongly on the value of w, with the 
(zi — 1.5, z a = 2.5) possessing changes of ~ 40% at w = 1 
as compared to w ~ 0, although, for the observed lensing 
systems, the same range in w produces a change of ~ 10% 
in the same quantity. Interestingly, for a fixed combination 
of Q m and f2 w , all the curves, irrespective of the redshift of 
the lens and source possess a minimum at the same value of 
w. The location of the minimum is tabulated in Table |l| 

At present, the general analytic solution to Equation [j] 
is quite complex (Giovi & Amendola 2001) and it is difficult 
to further analyze the minimum seen in Figures |lj For one 
case, where Q m = 0.0 and Q w = 1.0, however, analytic solu- 
tions for the relative angular diameter distance are straight 
forward. As the overall curvature is flat, the angular diam- 
eter distance between us and a distant source is 
D c (z) 



D A {z) = 



1 + 



where D c (z) is the comoving distance and is given by 
, dz' 

D.(z) = — 



HoJ E(z') 
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where E(z) = £L W (0)(1 + z)i {1+w) . Again, as the overall 
cosmology is flat, the angular diameter distance between 
the source and the lens is given by; 

D A {z 1 ,z 2 ) = —^—[D c (z 2 )-D c {z 1 )] (10) 

1 + 22 

With this, the relative time delays for differing values of w 
is seen to be: 



H [zi , z 2 ] oc 
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( 



where / = — i(l + 3w). This function possesses a minimum 
at w = — § , which is independent for the redshift pairs under 
consideration, as seen in Figure |l|. 

The preceding section has considered arbitrary combi- 
nations of cosmological parameters. Here, the impact of the 
choice of w on the estimation of H a in the favoured cosmo- 
logical model with Q m = 0.3 and flu, = 0.7 (top right panel 
in Figure |l|) is examined in more detail. Often, a cosmo- 
logical model where Q m — 1 is choosen when estimating H 
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Figure 2. The top panel presents the angular diameter distance of a source at z for an observer at z = 0, for the evolutionary w models 
discussed in the text. The thick lines correspond to universes with wo = — 1, while the thin line represents wo = —0.7. The solid line 
represents a non-evolving model with w± = 0, while the dashed lines evolve as w\ = 0.2 and the dot-dashed lines as w\ = 0.4. The lower 
panel presents the percentage difference between the evolving models and non-evolving model. 



from time delays, which is equivalent to considering w — in 
Figure |l[ An accelerating cosmology, as suggested from the 
cosmological supernovae programs, requires w < — i. It is 
interesting to note that in this cosmology, the minimum H a 
occurs at w = —0.15 (Table 1) the values of H derived at 
w — or w = — | will be very similar. Considering more neg- 
ative values of w, adopting classical cosmological constant, 
with w = — 1 results in H values that are ~ 5% different 
to the all matter case (w — 0). For several combinations 
of redshifts, however, (including the chosen fiducial points) 
more substatial differences of ~ 12 — 15% in the determined 
value of H - Considering w < — 1 results in more significant 
discrepancies, a point that will be returned to later. 



4 EVOLVING QUINTESSENCE 

One interesting aspect of quintessence cosmology is that, 
unlike classical Q m , cosmologies, it is possible that the 
equation of state w may vary over the cosmic history of the 
universe. Typically, a linear model where w(z) ~ wo + wjz 
has been adopted (e.g. Goliath et al. 2001). In the follow- 
ing analysis, it is assumed that O m = 0.3 and £l w = 0.7, 
consistent with the recent supernova experiments (Riess et 
al. 1999; Perlmutter et al. 1999) results. Four models for 
the evolution of w are adopted [taken from Linder (2001)], 



namely, (tuo.tui) = (-0.7, 0.2), (-0.7, 0.4), (-1.0, 0.2) and 
( — 1.0, 0.4). As the cosmology is flat, we follow the approach 
outlined in Equations H I and 0, with 



and 



f w (z) = exp 



+ w(z') 

1+2' 



dz 



(12) 



(13) 



Figures |2| and ^ present the angular diameter distance 
in these evolving models for an observer at z — and z = 1 
respectively. It is interesting to note that for an observer 
at redshifts greater than zero the fractional difference be- 
tween the angular diameter distances for an evolving and 
non-evolving quintessence component does not converge to 
zero as z s — > zi, rather, as seen in the lower box of Figure ^ 
they converge to a non-zero value. This seems somewhat 
counter-intuitive given the convergence of distances for an 
observer at redshift zero (Figure |^) . Employing Equation |l(] 
and setting Z2 = z\ + Az, the angular diameter distance 
becomes; 



D A (zi,zi + Az) 



Az dDJz) 



1 + Z! 



dz 



(14) 



where Az <C 1 + Z\. Therefore, the asymptotic ratio of the 
angular diameter distance in a non-evolving cosmology, with 
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Figure 3. As for Figure [3. but for an observer at z = 1.0. 
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Table 2. The ration of the time delay (or conversely Ho) for 
the cosmologies with linearly evolving quintessence. The first two 
columns are the lens and source redshifts, while the next four 
columns present the time delays for (wq,wi). Those with wq = 
— 1 arc normalized with respect to (—1.0,0.0), while those with 
wo = —0.7 are normalized with respect to (—0.7,0.0). 



a constant wq, and one where w(z) = Wo + zwi, as z-z — > Z\ 
simply becomes; 

-Pa (21, wo) _ j n m (l + zi) 3 + Mj + Zl f^o) 
where 

g(z) = exp {3(zwi + (1 + w a — wi) log[l + z})} . (16) 

With the values under consideration in Figure ^| this is pre- 
cisely the asymptote seen. 



Figure 4. The cosmological minimization function, oh / (H a ), 
versus the quintessence equation of state for the chosen fiducial 
cosmology and a sample of one hundred gravitational lens system. 
The dot-dash line represents the ideal case with no noise; as can be 
seen this is zero at w = — 1 . The black line assumes an uncertainty 
in the time-delay and lens modeling of 5%, the lighter curve of 
10% and the lightest 15%. A DC component has been subtracted 
from these latter curves. 
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As seen from Equation ^, the cosmological component 
of the measure of Hubble's constant is dependent upon a 
combination of angular diameter distances. Table |^ presents 
the relative values in the measured Hubble's constant for 
the linear evolution models described above, for the lens sys- 
tems in Section ^| The largest changes are seen in the first 
two rows which represent fiducial models; these possess the 
highest lens redshifts in the sample. For the observed lens 
system, it is clear that the evolving quintessence models do 
differ from non-evolving models, this difference is not sub- 
stantial, being typically less than a few percent for evolution 
models with ui\ = 0.4. 



5 MONTE CARLO SIMULATIONS 

How many lensed systems are required before firm limits can 
be made? To address this question a number of Monte Carlo 
simulations were undertaken with a population of lensed 
sources, assuming that for each a time delay could be de- 
termined to a specified accuracy. For real lensed systems, 
additional uncertainty is introduced from the lens modeling. 
Here, however, it is assumed that the lens modeling intro- 
duces no systematic uncertainty and that the random error 
is taken into account as an additional source of uncertainty 
in the time delay. 

A number of programs have examined the relative 
statistics of gravitational lensing (e.g. Kochanek 1993a; 
1993b; Williams 1998; Falco et al. 1998; Helbig et al. 1999). 
The complex analyses in these studies go beyond this pa- 
per and a more straight forward approach was adopted; this 
was to examine the distribution of source and lens redshift 
pairs in the CASTLES database (Munoz et al. 1998). This 
is seen to be a scatter plot of sources between redshifts 1 
and 4, with lenses between 0.2 and 1.1. Source (z s ) and lens 
(zi) redshifts where, therefore, selected randomly in these 
range, but ensuring that z s > 4z; — 2.2. This cut ensures 
that lens and sources are always reasonably separated in 
redshift space. 

In choosing a background cosmology, a fiducial model 
with £l m — 0.3 and £l w = 0.7 with a w = —1 and H = 1 
was employed. In an ideal universe, where there are no mea- 
surement errors and where gravitational lens models can 
be uniquely determined, the analysis of a sample of sys- 
tems may still produce a distribution of H values, as an 
incorrect choice of cosmological parameters, including the 
quintessence equation of state w, influences differently each 
system (e.g. Figure [j]). In this ideal universe, however, all 
one would need to do is to vary the cosmological parame- 
ters until all systems yielded the same value for the Hubble 
constant. In this way, not only would H be measured, but 
the underlying cosmology would be determined as well. 

In the real universe, however, the influence of measure- 
ment uncertainties needs to be considered. Instead of ob- 
taining a unique value of Hubble's constant, one could vary 
the cosmological parameters such that the dispersion in the 
resultant H 's are minimized. An example of this for a pop- 
ulation of one hundred gravitational lens systems is given 
in Figure ^; note that the function that is minimized is 
a h I {Ho} as it is the fractional dispersion of the results that 
is of interest. The dot-dashed line is the ideal universe case, 
where there are no sources of uncertainty, which possesses 
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Table 3. The 95% confidence intervals for w (upper values) and 
H (lower values) for the Monte Carlo samples described in the 
text. The sample sizes are given across the top of the table, and 
each row is for an assumed error in the determination of the time 
delay. Note that in calculating these confidence intervals, outlying 
points with w > have been neglected so that they represent the 
width of the distribution. 

a minimum value (of zero) at the chosen fiducial value. The 
solid curves represent varying values of noise in the determi- 
nation of the time delay, the black being 5% noise, grey 10% 
and light grey 15% noise. Clearly this function broadens as 
more noise is added, and the minimum is not necessarily at 
the fiducial value of w. 

Figure ^| presents the results of undertaking this pro- 
cedure for various samples of gravitational lens systems. 
The left panel presents the probability distribution for the 
quintessence equation of state w while the right hand panel 
is the corresponding distribution in Hubble's constants; note 
that all the distributions have been normalized to a peak 
value of one. The top row is for a sample of ten gravitational 
lens systems, akin to the situation today, followed by 50, 100 
and 500 systems. Ten thousand realizations were undertaken 
for each sample size. Each line on the plot corresponds to a 
different level of uncertainty in the gravitational lens time 
delay as outlined in the previous paragraph. Clearly today, 
where there are but a handful of gravitational lens systems 
for which we have determined the time delay, with overall 
uncertainties exceeding 5-10%, then the resultant Hubble's 
constant and w that can be derived from the sample are 
unlikely to accurately represent the underlying values. 

Table ^ presents the 95% confidence interval for the esti- 
mation of w and H for the various samples. Increasing the 
sample size greatly improves the situation, although even 
with 500 lenses and 15% noise, the values of w and H are 
not strongly constrained. One is led to conclude, therefore, 
that a large sample of lenses with very accurately determined 
time delays and lens models is required to significantly de- 
termine the the underlying cosmological parameters. Given 
the observational effort in such a task, other approaches to 
probing the cosmological equation of state are likely to prove 
more fruitful. Given this, the analysis was not extended to 
consider the smaller influences of quintessence evolution (see 
Section §). 

6 DO CURRENT LENS SYSTEMS SUGGEST 

W < -1? 

Can the current observations of time delays in gravitational 
lens systems tell us anything about the equation of state 
of the quintessence component? In recent years, dedicated 
monitoring of a number of lensed systems has provided ac- 
curate (~ 10 — 20%) time delay determinations (e.g. Fass- 
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Figure 5. The results of the Monte Carlo simulations of the determination of the quintessence equation of state and Hubblc's constant 
from the measurement of samples of gravitational lens time delays. The black line assumes an uncertainty in the time-delay and lens 
modeling of 5%, the lighter curve of 10% and the lightest 15% 



nacht et al. 1999; Koopmans et al. 2000). An examination 
of the Hubble's constant derived from such studies reveals 
that, typically, it is less than the value determined from lo- 
cal studies, even accounting for standard cosmological differ- 
ences (e.g. Impey et al. 1998; Koopmans et al. 2001; Winn et 
al. 2002). This very question was also recently addressed by 
Kochanck (2002) who suggests that this discrepancy is po- 
tentially due to galaxies possessing concentrated dark mat- 
ter halos with a constant mass-to-light ratios, at odds with 
expectations from cold dark matter structure models. Here, 
an alternative solution is considered. 

Examining the panel in Figure 1 corresponding to an 
fim = 0.3, fl w — 0.7 cosmology, it is apparent that choos- 
ing fi m = 1 (equivalent to setting w = 0) results in al- 
most the lowest possible determination of H D . Considering 
a cosmology with a classical cosmological constant (w = — 1) 



increases the determined value of H by ~ 5 — 20% (depen- 
dent upon the source and lens redshifts), but as noted above, 
the currently determined values still tend to lie below the 
72km/s/Mpc derived locally (Freedman et al. 2001). Assum- 
ing that the gravitational lens models are correct, one way 
to reach concordance between the two approaches is that 
w < — 1; such a conclusion is consistent with the recently de- 
rived limit of w < —0.85 from an analysis of a combination 
of cosmic microwave background, high redshift supernovae, 
cluster abundances and large scale structure data (Wang et 
al. 2000; Bean & Melchiorri 2001). While tantalizing, how- 
ever, it must be conceded that the current differences in the 
approaches are not statistically significant, especially given 
the relatively large uncertainty in the modeling of gravita- 
tional lens mass distributions (see, for example, the range 
of H values obtained from the modeling of PG 1115+080; 
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Kochanek 2002). Given a large sample of gravitational lens 
systems, however, the determination of a systematic differ- 
ence in the derived Hubble's constant could be made. 
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7 CONCLUSIONS 

This paper has investigated the role of a quintessence com- 
ponent on angular diameter distances, specifically their in- 
fluence on the determination of Hubble's constant from the 
measurement of a time delay in multiply imaged quasars. 

For flat universes, with an unevolving quintessence com- 
ponent, its seen that, for a gravitational lens system in which 
the time delay has been measured, the resultant Hubble's 
constant is dependent upon the value of the equation of 
state parameter w. Interestingly, the dependence of the de- 
termined value of the Hubble's constant as a function of w 
possesses a minimum which is independent on the lens and 
source redshift. 

Several models of evolving quintessence were also ex- 
amined, consisting of a linear evolution of the equation of 
state with redshift. The cosmologies resulted in significantly 
different forms of the angular diameter distance. Hence, our 
view of the cosmos would be different in the various cos- 
mologies. When considering the specific combination of an- 
gular diameter distances that constitute the cosmological 
contribution to the gravitational lensing determination of 
Hubble's constant, it is seen that the resulting variations 
between cosmologies is very small, a matter of only a few 
percent, relative to an unevolving case with the same present 
day constitution. 

A number of Monte Carlo simulations of the determi- 
nation of Hubble's constant and the quintessence equation 
of state, w, were undertaken to explore the efficacy of this 
approach. These revealed that the present situation with 
only a handful of lensed systems does not allow an accurate 
determination of the cosmic equation of state, and that at 
least an order of magnitude more lenses are truly required 
to provide a reasonably robust determination of the under- 
lying cosmology. The next generation of all-sky surveys are 
presently underway (e.g. Sloan Digital Sky Survey) or are 
being planned (e.g. VISTA, PRIME), and these datasets 
will greatly increase the number of lensed quasars available 
for monitoring studies. Cooray & Huterer (1999) estimate 
that ~ 2000 lensed quasars will be identified from the SDSS 
database alone, and a much larger number can be expected 
from the deeper VISTA and PRIME surveys. The number of 
these sources amenable for follow-up monitoring campaigns 
will naturally be much smaller, but one can confidently ex- 
pect a sample of several hundred systems to eventually be- 
come available. However, given the effort required to first 
find such systems, as well as monitor them to determine the 
time delays and the modeling procedure, it is likely that 
w will be first determined using one of the other various 
techniques currently being proposed. We conclude, there- 
fore, that gravitational lens time delays are likely to prove 
poor probes of the universal equation of state. 
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